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Abstract: The reaction pp — )• tbX is found to be rather rich in 
exhibiting several different types of CP asymmetries. The spin 
of the top quark plays an important role. Asymmetries are re- 
lated to form factors arising from radiative corrections of the tbW 
production vertex due to non-standard physics. As illustrations, 
effects are studied in two Higgs Doublet Models and in Super- 
symmetric Models; asymmetries up to a few % may be possible. 

The origin of CP violation remains a pressing issue in Particle Physics. 
The Standard Model (SM), with three generations of quarks, can accommo- 
date a CP violating phase, the Cabibbo-Kobayashi-Maskawa (CKM) phase 
[1]. However it is widely believed that this phase cannot account for baryoge- 
nesis [2]. Additional CP violating phases due to new physics are therefore a 
necessity. Besides, in extensions of the SM, new phase(s) appear rather read- 
ily. It is therefore quite unlikely that the CKM phase is the only CP violating 
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phase in nature. In particular, in top physics the SM causes neghgible CP 
violation effects [3] whereas, in sharp contrast, non-standard sources often 
give rise to appreciable effects [4, 5]. Searching for CP violation in top-quark 
production and/or decay is therefore one of the best ways to look for signals 
for new physics. 

In this paper we examine CP violation asymmetries in top quark produc- 
tion via the basic quark-level reaction: 

u + d^t + b, u + d^t + b (1) 

Indeed the reaction is very rich for CP violation studies as it exhibits many 
different types of asymmetries. Some of these involve the top spin. Therefore 
the ability to track the top spin through its decays becomes important and 
top decays have to be examined as well. 

It is easy to see that in the SM these effects are expected to be extremely 
small as they are severly suppressed by the GIM mechanism [3]. As illustra- 
tion of the possibilities with non-standard sources of CP violation we consider 
two extensions of the SM: a two Higgs Doublet Model (2IIDM) with Natural 
Flavor Conservation (NFC), often called type II model and a Supersymmet- 
ric Standard model (SSM). We find that CP asymmetries can be sizable, 
in some cases at the level of a few percent. Thus the asymmetries in tb 
production can be appreciably larger than those in tt pair production [6, 7] 
wherein they tend to be about a few tenths of percents. Since the number 
of events needed for observation scales as (asymmetry)"^ the enhanced CP 
violation effects in tbith) may make up for the reduced production rates for 
tb compared to tt. In fact larger asymmetries are not just gratifying but they 
can also be essential as detector systematics can be a serious limitation for 
asymmetries < 0.1%. 

Let us first discuss the asymmetries in the ud. (ud) subprocess. We con- 
sider four types of asymmetries that may be present. First is the CP violating 
asymmetry in the cross section: 

Ao = {(Tq - (Tq)/ {(Tq + CTq) (2) 

where (Tq and (Tq are the cross-sections for ud — )• tb and ud — )• tb respectively 
at s = [pt -\-piy. The CPT theorem of quantum field theory implies that the 
total cross-section for ud and du are identical. If a cross-section asymmetry 



2 



Aq in the tb final state is present, then to maintain the balance of total cross 
sections, another mode must have a compensating asymmetry. 

The spin of the top allows three additional types of CP violating polar- 
ization asymmetries. To define these let us introduce the co-ordinate system 
in the top-quark rest frame where the unit vectors are given by oc — Pfe, 
Cy oc Pu X Pf, and = Cy X e^. Here Pj, and P„ are the 3-momenta of the 
b quark and the initial u-quark in that frame. We denote the longitudinal 
polarization or helicity asymmetry as 

A{z) = [Nn -Nl-Nl + Nn)/{Nn + Nl + Nl + Nn) (3) 

where Nl is the number of left-handed top quarks produced in ud — )• tb and 
Nr is the number of right-handed t produced in ud — )• tb etc. In the frame 
introduced above therefore right handed tops have spin up along the 2;-axis 
and left handed ones spin down. We further define the CP violating spin 
asymmetries in the x and y directions as follows: 

A{x) = {N,+ -N,_+N,+ -N,_)/{N,+ + N,_ + N,+ + N,_) (4) 

A{y) = {Ny+-Ny_-Ny++Ny_)/{Ny++Ny_+Ny++Ny_) (5) 

where Nj^[Nj^) represent the number of t[t) with spin up with respect to j 
axis, for j = x, y. 

While all these four asymmetries are manifestly CP-violating, Aq, A[z) 
and A[x) are even under naive time reversal (Tjv) whereas A[y) is Tjv-odd. 
So the first three require a complex Feynman amplitude whereas A[y) needs 
a real amplitude. Of course, all four do need a CP-violating phase in the 
underlying theory. 

In the limit of massless u and d quarks the CP violating contribution to 
the Wtb vertex may be represented by the effective interaction: 

C = i2-^gw W+ t[F-f^ + im;^Ga"^qa\ Lb 

- i2-^gw W- b R [F-f^ + im;^Ga"^qa\ t (6) 

where L = [1 — 75)/2, R = 75)/2, we have taken the b-quark to be mass- 
less and consider only the left handed component which will interfere with 
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the standard model (tree-level) contribution. The asymmetries of interest 
are therefore expressible in terms of F and G alone. We will denote the real 
and imaginary parts hy F = Fr + iFj and G = Gr + iGj. Note that the 
terms proportional to Fr and Gr are hermitian while the terms proportional 
to Fi and Gj are not and therefore proportional to final state interaction 
effects. 

We can now derive expressions for the asymmetries in terms of the form 
factors, for example: 



^0 = -2F, + -^Gn A{z) = 2\-^F, - -^G: (7) 
z + x z + x z + x 

where x = m^/s. Notice that the three quantities {Aq, A[z)^A[x)} are 
linear combinations of the two form factors {F/, G/}, so in particular, one 
can show that 

A{x) = -3Trx-2 ((2 + x)Ao + (2 - x)A{z)) /32 (8) 

The dependence of Aq on Fj and Gi provides a clue as to how the balance of 
total cross-section required by CPT is achieved. In order for these imaginary 
parts to exist in perturbation theory, there must be a contribution from a 
loop graph which has an intermediate state J that can be kinematically on 
shell. J is therefore another component of the cross-section, and in fact it is 
the cross-section asymmetry in J that compensates for the asymmetry in tb. 

The asymmetry on the other hand, is proportional to Gr. Since 

Gr may be related to Gi through a dispersion relation (as a function of s) 
we can obtain A[y) in terms of Aq and A[z): 

A{y)[s 




■m - ^)Mm + {2c + x)A{msM^ (9) 

Note that the integrand is if is below the threshold to produce an imag- 
inary part. In deriving this relation, we assume that some GIM like cance- 
lation applies in the underlying model in which case all of the asymmetries 
vanish in the limit of large s. 
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Fig. la shows the SM tree-level production process. The necessary ab- 
sorptive parts require radiative corrections, involving a CP violating phase, 
at least to one-loop order. Fig. lb shows the only graph relevant to a type-II 
2HDM with a CP phase residing in neutral Higgs exchanges. Fig. Ic shows 
example of a one-loop graph that pertains to a SSM which can involve new 
CP violating phase(s) as well as the needed absorptive parts. 

As is well known, in a 2HDM with NFC, CP nonconservation emanates 
from soft symmetry breaking complex parameters in the Higgs potential [8, 
9]. These induce mixing between real and imaginary parts of the Higgs 
fields in their mass matrix. Consequently the mass eigenstates do not have 
a definite CP property. Therefore, an important feature of the 2HDM is 
that CP-violation may result from the neutral Higgs sector even when there 
is none in the charged Higgs sector. This CP violating phase from neutral 
Higgs exchanges is much more difficult to look for compared to that from the 
charged Higgs exchanges. The top quark can play a special role with regard 
to the neutral Higgs CP as due to its large mass its coupling with the Higgs 
are significantly enhanced compared to all the other quarks. 

A distinctive manifestation of such CP violation is that the neutral Higgs 
mass eigenstates couple to fermions with both scalar and pseudoscalar cou- 
plings. The relevant Feynman rules for calculating the asymmetries of inter- 
est can be extracted from the following part of the Lagrangian, 

£^0 = i?°/(aj, + z6/,75)/ + cwjmwH'lg,,W''W'' (10) 

This involves the ffH^ and the WW couplings, where j = 1,2,3 for the 
three neutral spin fields. The coupling constants, ajj, hj^ and are 
functions of tan /3, which is the ratio between the two vacuum expectation 
values (VFVs) in this model, i.e. tan /3 = ^2/^^!, and of the three mixing 
angles ai...3 which diagonalize the 3x3 Higgs mass matrix. 

For simplicity we assume that two of the three neutral Higgs particles 
are much heavier compared to the third one. The effects we are seeking are 
therefore likely to be dominated by the lightest neutral Higgs. We thus omit 
the index j in Fqn. 10 and denote the couplings of the lightest Higgs with 
the top and the W a,s at, bf and cw- 

From Fig. 1 we see that the imaginary part of the loop is provided by 
the WH intermediate state and hence the cross-section asymmetry Aq is 
compensated by an asymmetry in ud — )• W'^ H versus ud — )• W~ H . Clearly 
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this imaginary part can only exist above the WH threshold at s = {niw + 
muY . Thus below this Aq, A{z) and A{x) will be identically though A{y) 
need not be since it depends only on virtual effects. 

Using the Lagrangian (fO) the CP asymmetries, Aq and A{z), resulting 
from the interference of Fig. la and lb can be readily calculated: 

Ao = -hS^p^Ul -3y- z)<P - 2(1 - y){x + xy - xz - iy)r\ (11) 



A{z) = Ji£ll!^((l + 3.-7,-. + 3., + ..)<^-2r 
lD7rmt(l — X) [ 

-\-[3x — 4j/)(l — z -\- xz) + x[l — x)y[y — z) 



{x-2yf 
(121 



where x = m'^/s, y = m^/s, z = mfj/s, cj) = \/\ -\- y^ -\- z^ — 2y — 2z — 2yz, 
Rq = x/[y[2 + x)[l — x)) and 

T = (1 — a;)~"^tanh~"^ [(1 — x)(f)/[l -\-x — y — z-\-xz — xy)] (13) 

As should be expected from Eqn. (10) all the CP asymmetries have to 
be proportional to the product htCw ■ We choose the angles in the Higgs 
mixing matrix as cti = a2 = 7r/2 and cts = as it tends to give maximal 
effects[10, 11]. Then one can show that bfCwmw — .2mt cos /3 cot /3 and the 
asymmetries can now be expressed as a function of tan /3 and niH only. 

We present our numerical results for tan /3 = 0.3 [12]. Numbers for other 
values of tan /3 can then be readily obtained using the relation given in the 
preceding paragraph. Fig. 2 a and b show the asymmetries as a function of 
s for niff = 100 and 400 GeV. The cross section asymmetry, Aq, and helicity 
asymmetry, A[z) are typically a few percent. The polarization asymmetry 
A[x) likewise is a few percent. 

Since the real part of the graph in Fig. lb does not need a physical 
threshold, it may receive contributions from Higgs bosons of arbitrary mass. 
In the limit of degenerate Higgs masses, CP violating effects should vanish. 
Hence it may not be a legitimate approximation to ignore the more massive 
Higgs bosons in this loop. We will assume therefore that the other Higgs 
bosons of the theory have a mass and, for the purposes of numerical 
estimates we will take = ITeV. The effect of these heavy Higgs bosons 
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on the dispersion Eqn. (9) is to replace {Aq, by {Aq, which are 

defined by subtracting the efi'ects of the heavy Higgs. So, for instance, Aq 
gets replaced by A'q[s] = Ao[s,m//] — Ao[s,m^] and likewise for the other 
asymmetries. Note that if s < (m^ + niwY then A'q = Aq; A'(i, x) = A{z, x) 
since s is below the W , H' threshold. In Fig. 2a and 2b we also show the 
value of A{y). There tends to be a cusp at the HW threshold where the peak 
value of A{y) is about one percent. 

To search for the effects of these three types of spin asymmetries that 
occur at the production vertex, decays of the top will obviously need to 
be examined [13]. In particular when considering A{y) one must keep in 
mind that it is dependent on the real part of the loop amplitude of Fig. lb. 
One complication that this could lead to is that a similar asymmetry may 
also enter into the decay t — )• hW when similar radiative corrections to that 
vertex are also included [15]. This is not a concern in the case of the other 
observables since if we assume that the Higgs is above the threshold, i.e. 
{niw + niH) > rrit, the necessary condition that there be an imaginary part 
in the decay amplitude is not satisfied. 

As it turns out, the observed value of A{y) is not affected by CP vi- 
olation in the decay process. The key point is that the measurement of 
A{y) through the decay chain u[pu) d[pd) — > b[pbi) t[pt) followed by t[pt) — > 
b{pb2) e~'~(pe) i^iPiy) is equivalent to measurement of the term proportional 
to e{pe,pd,pt,Phi)- On the other hand, CP violation arising from the decay 
process is proportional to e(pe, Pd, Pt, ^62)- It is easy to see that an observable 
related to the first of these will be insensitive to the second[IO]. 

These quark level asymmetries can be readily converted to the hadron (i.e. 
pp) level by folding in the structure functions in the standard manner [16]. 
The results are shown in Fig. 3 where for each of the asymmetries we apply 
a cut of s > {niH + niwY ■ At the Tevatron [E = 2 TeV) the expected 
number of events are 900-3000 with an integrated luminosity 3/6~^-I0/6~^ 
respectively[I7]. If the collider energy gets upgraded to 4 TeV and/or there 
are additional luminosity upgrades as have often been discussed, then the 
number of events can go up by another factor of about 2 to I0[I7]. Thus the 
asymmetries, in the range of a few percent, resulting from some extensions of 
the SM may well become within the reach of experiment provided, of course, 
that the signal for these single top events could be extracted from possible 
backgrounds[I8]. 

Another extension of the standard model which can produce these kind 
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of asymmetries is SSM. There are a number of possible graphs which could 
contribute[10, 11]; here we will consider only the gluino exchange diagram 
given in Fig. Ic. In this case CP violation arises through the mixing matrix 
between the fermion and the scalar states, in general a 6 X 6 matrix. For sim- 
plicity, let us consider a scenario, motivated by supergravity models, where 
all the squarks are degenerate with a mass rhq except for the super partner of 
the top quark, the stop. Furthermore, the two stop states mix with the left 
and right parts of the top quark with a general 2x2 unitary mixing matrix 
A! . The cross-section asymmetry therefore is given by: 

Ao = 2aJm{XiiXl2) [g{rhti, rhb, rhg, s) - g{mt2, rhb, rhg, s)] (14) 

where rhti^2 are the masses of the two stops. The function g is given by: 

g = (4/3)7^^(1 - x)-^{2 + x)-V(Q"^tanh-^g - 1) (15) 

where x = m^/s^ y = rh'^/s^ z = rh^/s and w = rh'^/s while (f> has the same 
definitions as previously in terms of y and z and Q = [1 — x)(f> / [1 -\- 2w — x — 
y — z -\- xz — xy) 

In this case the helicity structure of the model is such that the form 
factor F = G. Thus A{z) = Aq from which A{x) and A{y) follow through 
Fqn. (8,9). In Fig. 2 we also show these asymmetries due to the SSM 
for mti = 100 GeV, mt2 = 500 GeV, = 100 GeV and = 100 GeV. 
We have also assumed that the the quantity Im{XiiX22) = 1/2 which is 
its maximum value. We can see that in this case the asymmetries are less 
than 1%. The small size of these asymmetries is, in part, due to the fact 
that the intermediate state (see Fig. Ic) consists of two scalars that must 
be in a P-wave giving rise to an additional threshold suppression factor. 
However, in SSM many other types of loop corrections (e.g. box graphs) can 
also contribute giving rise to asymmetries on the order of several percents. 
[10, 11]. 

We close with a few remarks in brief. 

1. It is important to note that from the point of view of experimental 
detection these four asymmetries are independent. Thus, the sensitivity 
of a given detector to observing the combined CP violation effects may 
be appreciably better than that for any one asymmetry [10, 11]. 
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2. We have focussed here on a pp machine (i.e. the Tevatron) as the self- 
conjugate nature of the intial state is rather important for CP studies. 
At the LHC ( i.e. a pp machine), although the event rate is high, such 
CP studies are quite difficult. Note, for instance, that the cross-sections 
for pp — )• tbX and to tbX are expected to be different at the LHC even 
if CP was strictly conserved. 

3. We recall that the l^-glue fusion subprocess, W'^ -\- gluon — )• t -|- 6, 
also contributes to the same final state[f7]. While it will be useful to 
include its contribution to the asymmetries in a future study, for now 
we note that, at least in the 2HDM, CP violating radiative corrections, 
to one loop order, to Wg fusion do not yield absorptive parts (in the 
nib = limit). 

We will address to some of these issues in greater detail in future work 
[10]. 

This research was supported in part by the US-Israel Binational 
Science Foundation and in part by USDOE contracts DC-AC05- 
84ER40I50 (CEBAF) and DE-AC-76CH000I6 (BNL). 
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Figure Captions 



Figure 1: Feynman diagrams for contributions to ud — )• tb. la) the stan- 
dard model process, lb) one-loop graph in the two Higgs doublet models 
(2HDM), and Ic) an example of one-loop graph that could occur in the 
SSM. 

Figure 2: The magnitudes of the quark-level asymmetries Aq (solid); A[z) 
(dashed); A{x) (dotted) and A{y) (dot-dashed) as a function of \/j in the 
(2HDM) with tan /3 = 0.3. Fig. 2a: niH = 100 GeV and Fig. 2b: niH = 400 
GeV. Note that A[y) is computed keeping fixed the masses of the two heavier 
neutral i7°'s to 1 TeV. Also shown with the lower solid line is the asymmetry 
Aq in the SSM described in the text for parameters rhti = lOOGeV, mt2 = 
bOOGeV, rhg = lOOCeV, rhg = lOOGeV and ImiXuX^^) = 1/2. 

Figure 3: The corresponding asymmetries in the pp cm., for E = 2 TeV, 
frame as a function of niH for the 2HDM and as a function of Amt(= mt2 — 
rhti) for the SSM case. See also caption for Fig. 2. 
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Figure 2b 
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